AN 84-/iG MAGNETIC FIELD IN A GALAXY AT Z=0.692 
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The magnetic field pervading our Galaxy is a crucial constituent of the 
interstellar medium: it mediates the dynamics of interstellar clouds, the energy 
density of cosmic rays, and the formation of stars l . The field associated with 
ionized interstellar gas has been determined through observations of pulsars in 
our Galaxy. Radio-frequency measurements of pulse dispersion and the rotation 
of the plane of linear polarization, i.e., Faraday rotation, yield an average 
value B « 3 /iG (ref. 2). The possible detection of Faraday rotation of linearly 
polarized photons emitted by high-redshift quasars 3 suggests similar magnetic 
fields are present in foreground galaxies with redshifts z > 1. As Faraday 
rotation alone, however, determines neither the magnitude nor the redshift of 
the magnetic field, the strength of galactic magnetic fields at redshifts z > 
remains uncertain. Here we report a measurement of a magnetic field of B m 
84 /iG in a galaxy at z =0.692, using the same Zeeman-splitting technique that 
revealed an average value of B = 6 /iG in the neutral interstellar gas of our 
Galaxy 4 . This is unexpected, as the leading theory of magnetic field generation, 
the mean-field dynamo model, predicts large-scale magnetic fields to be weaker 
in the past rather than stronger. 5 

We detected Zeeman splitting of the z=0.692, 21 cm absorption line in the direction 
of the quasar 3C 286 (refs. 6,7) using the 100-m Robert C. Byrd Green Bank Telescope 
(GBT) of the National Radio Astronomy Observatory. The absorption arises in a 
damped Lya system (henceforth denoted DLA-3C286) that is drawn from a population 
of neutral gas layers widely thought to be the progenitors of modern galaxies 8 . The 
radio data for DLA-3C286 are summarized in Fig. [TJ which shows the line-depth spectra 
constructed from the measurable quantities used to describe polarized radiation, that 
is, the Stokes parameters. We show the line-depth spectra constructed from the I(y) 
and V(y) Stokes parameters (where v denotes frequency) near the 839.4 MHz frequency 
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centroid of the redshifted 21 cm absorption line. Fig. [T^l shows the line-depth spectrum 
constructed from 1(f). A Gaussian fit to the absorption line in Fig. [Th^ yields a redshift 
^=0.6921526±0.0000008, central optical depth of r =0.095±0.006, and a velocity dispersion 
<r„=3.75±0.20 km s _1 , which are in good agreement with previous results 6,7 . 

In Fig. [Tb we plot the line-depth spectrum constructed from V(v), which shows the 
classic 'S curve' pattern expected for Zeeman splitting. From our least squares fit to the 
data, we find -Bi os =83.9±8.8 /iG, where B\ os is the magnetic field component projected 
along the line of sight (we note that the direction of B\ os is unknown since the instrumental 
sense of circular polarization was not calibrated). This magnetic field differs in two respects 
from the magnetic fields obtained from Zeeman splitting arising in interstellar clouds in the 
Galaxy. First, the field strength corresponds to the line-of-sight component of the mean 
field < -Bios > averaged over transverse dimensions exceeding 200 pc, as very-long-baseline 
interferometry observations of the 21-cm absorption line show that the gas layer must must 
extend across more than 0.03 " to explain the difference between the velocity centroids of 
the fringe amplitude and phase-shift spectra 9 (although the data are consistent with a 
-B-field coherence length less than 200 pc, the resulting gradient in magnetic pressure would 
produce velocity differences exceeding the shift of ~ 3 km s" 1 across 200 pc detected by 
very-long-baseline interferometry). By contrast, the transverse dimensions of radio beams 
subtended at neutral interstellar clouds in the Galaxy are typically less than 1 pc. Second, 
this field is at least an order of magnitude stronger than the 6-/zG average of B fields 
inferred from Zeeman splitting for such clouds 4 . 

We obtained further information about conditions in the absorbing gas in DLA-3C286 
from accurate optical spectra acquired with the HIRES Echelle spectrograph on the Keck I 
10 m telescope. Fig. [2] shows velocity profiles for several resonance absorption lines arising 
from dominant low-ionization states of abundant elements. The results of our least squares 
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fit of Voigt profiles to the data are shown in Table 1, where the optical redshift is displaced 
+3.8 ±0.2km s _1 from the 21 cm redshift. This solution also yields ionic column densities 
from which we derived the logarithmic metal abundances with respect to solar abundances, 
[M/H], and dust-to-gas ratios with respect to the Galactic interstellar medium, [D/G]. 
These are among the lowest values of [M/H] and [D/G] deduced for damped Lya systems 
at z=0.7 (refs. 10, 11). The low metallicity indicates a history of low star formation rates 
(SFR). Because the intensity of far ultra-violet radiation emitted by young massive stars is 
proportional to the concurrent SFR per unit area S SFR , low values of Ssfr should result 
in low grain photoelectric heating rates per H atom, r pe (ref. 11). This is consistent with 
the low upper limit, r pc < 1CT 27 ' 4 ergs s _1 per hydrogen atom, obtained by combining the 
assumption of thermal balance with the absence of C II* absorption (that is, absorption 
from C II in the excited 2 P^/2 fine structure state) at a wavelength of 1335.7 A in the 
previous low-resolution Hubble Space Telescope (HST) spectra of 3C 286 (ref. 12), and 
indicates Ssfr < 10~ 29 M Q yr -1 kpc~ 2 (95 % confidence level), which is less than the 
solar-neighborhood value of 10~ 2 - 4 M yr _1 kpc~ 2 (ref. 13). 

As a result, we have detected an unusually strong magnetic field at z — 0.692 with a 
coherence length that probably exceeds 200 pc in neutral gas that is quiescent, metal-poor, 
nearly dust-free, and presents little evidence for star formation. To model this configuration, 
we first consider the magnetostatic equilibrium of a plane-parallel sheet with in-plane 
magnetic field -B p i an e orthogonal to the vertical gravitational field exerted by gas with 
perpendicular mass surface density, S. In magnetostatic equilibrium the total midplane 
pressure, (i? 2 lane /87r) + per 2 , equals the 'weight' of the gas, 7rGS 2 /2, where p is the mass 
volume density of the gas and G is the gravitational constant. However, because the 
pressure-to-weight ratio exceeds 715 in DLA-3C286, the magnetized gas cannot be confined 
by its self-gravity. Therefore, self-consistent magnetostatic configurations are ruled out 
unless the contribution of stars to £ exceeds ~ 350 M Q pc~ 2 . Although this is larger than 
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the 50 M pc~ 2 surface density perpendicular to the solar neighborhood, such surface 
densities are common in the central regions of galaxies. In fact high surface densities of stars 
probably confine the highly magnetized gas in the nuclear rings of barred spirals. These 
exhibit total field strengths of ~ 100 /xG, inferred by assuming equipartition of magnetic 
and cosmic-ray energy densities 1 . However, because the rings are associated with regions of 
active star formation, high molecular content, and high dust content, they are unlikely sites 
of the magnetic field detected in DLA-3C286. 

On the other hand the absorption site might consist of highly magnetized gas confined 
by the gravity exerted by a disk of old stars. The H I disks found at the centers of 
early-type SO and elliptical galaxies 14 are possible prototypes. Support for this idea stems 
from a high-resolution image obtained with the Hubble Space Telescope: a Wide Field and 
Planetary Camera 2 (WFPC2) J-band image, from which the quasar has been subtracted, 
reveals residual emission spread over angular scales of ~ 1 " (ref. 15). The asymmetry of 
the light distribution with respect to the point-source quasar suggests that some of the 
light is emitted by a foreground galaxy with a brightness centroid displaced less than 0.5 " 
from the quasar. The location of diffuse emission in the direction of an amorphous object 
detected 2.5 " from the quasar in ground-based imaging 16 further suggests that the diffuse 
emission comes from central regions of the amorphous object. A recent reanalysis of the 
WFPC2 image shows the amorphous object to be a filament resembling a spiral arm or 
tidal tail (H.-W. Chen, personal communication), that is, the outer appendage of a galaxy 
centered within a few kpc of the quasar sightline. 

However, the magnetic field detected in DLA-3C286 may not be confined by gravity 
in an equilibrium configuration. Rather, the detected field may be enhanced by a shock 
(F.H.Shu, personal communication). Assuming a typical value of -B p i anc ~ 5 /iG for the 
equilibrium field of the preshock gas, we find that a shock-front velocity of ~ 250 km s _1 
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will result in a post-shock field strength of p=* 100 /xG in the limit of flux freezing in a 
radiative shock with post-shock density of ps 10 cm -3 . This scenario seems plausible 
because 250 km s -1 is reasonable for the impact velocity generated by the merger between 
the gaseous disks of two late-type galaxies, and the WFPC2 image is consistent with the 
presence of two foreground galaxies. But the second disk would create another set of 
absorption lines displaced > 250 km s -1 from the redshift of DLA-3C286, which is the only 
redshift observed. By contrast, the merger between a gaseous disk and an elliptical galaxy 
could result in only one damped Lja system redshift, as a significant fraction of ellipticals 
do not contain H I disks 14 . In this case a shock front moving in the plane of the disk galaxy 
would be generated by the gravitational impulse induced by the elliptical moving normal to 
the plane. Preliminary estimates indicate that an elliptical with a modest mass, M=2xlO n 
M Q , and impact velocity of ps 300 km s _1 would produce a cylindrical shock of sufficient 
strength to boost an initial field with -E> p i anc ~ 10 /xG to ps 100 /xG. 

Let us examine these scenarios more closely. The quiescent velocity field of the gas fits 
in naturally with the 'magnetostatic equilibrium' scenario, because the low value of Esfr 
suggests a low rate of energy injection into the gas by supernovae 17 , which could result 
in a velocity dispersion of a v 4 km s -1 . Moreover, the weak radio jets associated with 
early-type galaxies containing central H I disks are natural sources of B fields for these 
disks. However, 21 cm absorption measurements of such disks in nearby galaxies reveal the 
presence of absorption-line widths far broader than the narrow-line width of DLA-3C286 
(ref. 18). Also, it is unclear whether or not the high surface density of old stars required to 
confine the B fields are present in these disks, and whether or not the build-up of -Bpianc to 
100 /xG is possible in the 4 to 5 Gyr age of the disk. In the 'merger scenario' the dynamo 
need only build up to ps 10 /xG in the same time interval, but it is then necessary to explain 
why the post-shock velocity field averaged over length scales of 200 pc is so quiescent. 
Furthermore, the probability p for detecting ~ 100-/xG magnetic fields in a random sample 
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of 21 cm absorbers is small. Our estimates, based on the merger fraction of galaxies with z 
~ 1 (ref. 19) and on the duration time for magnetic field enhancement, suggest that p ~ 
0.005 to 0.03: either we were lucky, or some characteristic of DLA-3C286, such as narrow 
line width, is a signature of strong magnetic fields. 

Therefore, it is premature to decide among these and other possible models to explain 
the presence of the 84 /zG magnetic field in DLA-3C286. However, our data support 
the inference from recent tenetative evidence for Faraday rotation in high-z quasars 20 
that magnetic fields are generic features of galaxies at high redshifts, which potentially 
have a more important role in galaxy formation and evolution 21 than hitherto realized. 
Specifically, the highly magnetized gas that we have detected could suppress gravitational 
collapse and hence may be a reason for the low in situ star formation rates of high- ,2 
DLAs 22 . We plan to test this hypothesis by using the GBT to search for Zeeman splitting 
in high-redshift DLAs exhibiting 21 cm absorption. 
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Redshift ^=0.69217485±0.00000058 



Velocity Dispersion a v 


=3.08±0.13 km s" 1 


Ion X 


logioJV(X) 


[X/H] 




cm~ 2 




H I 


21.25±0.02 




Fe II 


15.09±0.01 


-1.66±0.02 


Cr II 


13.44±0.01 


-1.48±0.02 


Zn II 


12.53±0.03 


-1.39±0.03 


Si II 


> 15.48 


> — 1.31 



Table 1: Physical Parameters of DLA-3C286 Inferred from Optical Absorption 
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Fig. 1.- Line-depth spectra of Stokes-parameters. Data acquired in 12.6 hrs. of on-source integration with the GBT radio antenna. 
Because the GBT feeds detect only orthogonal linearly polarized signals, while Zccman splitting requires measuring circular polarization to construct 
Stokes V(f), we generated V(i^) by cross-correlation techniques 23 . The velocity v — km s~ 1 corresponds to z = 0. 6921526. a, Line-depth function 
Dj{v) = [I(v)— I c (v)]/I c (v), where I(v)= sq + sqq, with sq the power measured in linear polarization position angle 9, corresponds to the total 
intensity spectrum, and I c {v) is a model fit to the continuum. Dj(v) =exp[— t(l/)] — 1 where r(i^) = [T(f)Q + t"(i^)oq]/2 is the average optical 

depth in the two orthogonal states of linear polarization 4 . b, Line- depth function Dy (v) = V {v) / 1 c {v) , where V(v) =s^Qp — S LCP ^ s difference in 
power between the right-hand and left-hand circularly polarized (respectively RCP and LCP) signals. Here Dy{u) — — [ry (i/)/2]cxp[ — r{u)] , where 
Ty(f)=T R cp(£/) — r L cp(i') <<1 (ref. 4) is the difference between the optical depths of RCP and LCP photons. For Zeeman splitting of the 21 cm 
line, the degeneracy of the F — to F — 1 hypcrfinc transition is removed since the mp — —1, 0, +1 states differ in energy. This results in a small 
frequency difference between absorbed LCP photons (mjr — —1) and RCP photons (mp — +1). is crucial for detecting Zccman splitting 

because the orthogonal, circularly polarized states of the photon arc cigenstates of the spin angular momentum operator with eigenvalues ± h, 
that is, angular momenta directed along or opposite to the direction of photon propagation. 24 When B\ os = B, transitions between the hypcrfinc 
F — and F— 1 states occur exclusively through absorption of LCP or RCP photons through excitation of the mp = — 1 and mp — + 1 hypcrfinc 
states respectively. Because V (v) is the difference in the RCP and LCP intensities, the resulting V(v) line profile is the difference between two 
Gaussian absorption profiles with frequency centroids shifted by — 2.8(.Bi os //iG)(l + z)" 1 Hz (where Bi os is measured in microgaussj. The 

'S curve' is due to the sign flip in RCP-minus-LCP intensity difference as v passes through line centre. 
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Rdative Velocity (kms ) 



riff- — 

HIRES velocity profiles for dominant low-ionization states of abundant elements in the 21 cm absorber towards 3C 
286. Spectral resolution Au— 7.0 km s - " 1 and the average signal-to- noise ratio per 2.1 km s — 1 pixel is about 30:1. The bold dashed vertical 
line denotes the velocity ccntroid of single-dish 21 cm absorption feature and the faint dashed vertical line denotes the velocity centroid of the 
resonance lines shown in the figure. Our least squares fit of Voigt profiles (red) to the data (black) yields ionic column densities as well as the 
rcdshift ccntroid and velocity dispersion shown in Table 1 (lower and upper green horizontal lines refer to zero and unit normalized fluxes). Because 
refractory elements such as Fe and Cr can be depleted onto dust grains 2 ^ , we used the volatile elements Si and Zn to derive a logarithmic metal 
abundance with respect to solar abundances of [M/H] — — 1.30. The depiction ratios [Fc/Si] and [Cr/Zn] were then used to derive a conservative 
upper limit on the logarithmic dust-to-gas ratio relative to Galactic values of [D/G] < — 1.8. 



